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The masses of the pseudo-Nambu-Goldstone bosons in the color superconducting phase of dense 
QCD with two light flavors are estimated by making use of the Cornwall- Jackiw-Tomboulis effective 
action. Parametrically, the masses of the doublet and antidoublet bosons are suppressed by a power 
of the coupling constant as compared to the value of the superconducting gap. This is qualitatively 
different from the mass expression for the singlet pseudo-Nambu-Goldstone boson, resulting from 
non-perturbative effects. It is argued that the (anti-) doublet pseudo-Nambu-Goldstone bosons form 
colorless [with respect to the unbroken SU(2) C ] charmonium-like bound states. The corresponding 
binding energy is also estimated. 

11. 15. Ex, 12.38.Aw, 12.38.-t, 26.60,+c 

> 

^ . I. INTRODUCTION 

-xf . 

Although there is no reliable observational signature yet, from theoretical point of view, it is quite reasonable to 
assume that the cores of compact stars are made of color superconducting quark matter |lj . If we take this assumption 
seriously, it becomes quite important to study the properties of the possible color superconducting phases in full detail 
' (for an up to date review on color superconductivity see Refs. §,|J). In this paper, we continue our study of the 
Oh, pseudo-NG bosons, related to the approximate axial color symmetry, in the S2C phase of cold dense QCD. 

Let us recall that the axial color transformation is not a symmetry of the QCD action. However, its explicit 
[ breaking is a weak effect at sufficiently high quark densities where the coupling constant a s (/i) is small (p is a 
i chemical potential). This was our main argument in Refs. jl]|j| suggesting the existence of five (rather than one ||) 
light pseudo-NG bosons in the S2C phase of cold dense QCD. No mass estimates for these pseudo-NG bosons were 
given in Refs. @,|^. In this paper, we fill in the gap by developing the formalism and calculating the masses. 

we 
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This paper is organized as follows. In the next section, we briefly introduce our model and notation. In Sec. [II 



describe our method, based on the Cornwall-Jackiw-Tomboulis (CJT) effective action, for calculating mass estimates 



of the pseudo-NG bosons. Then, in Sec. IV, the leading order diagram is approximately calculated, using analytical 
methods. The fate of the colored pseudo-NG bosons in the doublet and antidoublet channels is discussed in Sec. M. In 
Sec. [v], we give our conclusions. In Appendices ^ and [B| we present the general expression for the gluon polarization 
tensor and the calculation of the integrals that appear in its definition. In Appendix C the problem of the gauge 
invariance in the loop expansion of the CJT effective action is discussed. 



II. THE MODEL AND NOTATION 



Here we consider cold dense QCD with two light quark flavors (u and d) in the fundamental representation of the 
SU(3) C color gauge group. In order to keep the analytical calculation under control, we assume that the chemical 
potential /i is much larger than Aqcd- Of course, when we talk about the compact stars, this is a far stretched 
assumption. Therefore, while extending our analytical results to the realistic densities existing, for example, at the 
cores of compact stars, one should be very careful. While most of the qualitative results may survive without being 
affected, most of the quantitative estimates would probably be valid only up to an order of magnitude. From the 
viewpoint of a theorist, it is still most interesting to study the predictions of the microscopic theory, i.e. QCD in the 
problem at hand. The price for such a luxury is the necessity to work at asymptotically large densities. 
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Instead of working with the standard four component Dirac spinors, in our analysis below, it is convenient to 
introduce the following eight component Majorana spinors: 



4% = Ci> 



D ) 



(1) 



where ipo is a Dirac spinor and C is a charge conjugation matrix, defined by C 1 7 /i C = — 7J and C = —C T . In this 
notation, the inverse fermion propagator in the color superconducting phase reads ||^-[Tl|] 

(Gwr 1 = - 



(p + /i)7°+/ A 

A (po-A*)7°+^ 
= _ { ( 7° [(Po ~ e~ )A+ + (po + e+)A-] A 

V A 7 [{ PQ - e+)A+ + {po + e-)A~ 

where A = 7°A t 7°, A% = ~f 5 e ij e ab3 [A~A" + A+A+] and the "on-shcll" projectors A^ ) ±) are 



p 2 ^ \p\ 



, where a = 7°7, 



(2) 



(3) 



(note that A* are complex valued gap functions). Color and flavor indices are denoted by small latin letters from the 
beginning and the middle of the alphabet, respectively. By definition, e~ = \p\ ± /1 and f = —p - 7. The effects of the 
quark wave function renormalization are neglected here. 

Now, after inverting the expression in Eq. (0), we arrive at the following propagator: 



G( P ) 



Ru R12 
R21 R22 



where 



Rn(p) = J°li 
R 22 (p) =7°^i 
Ruip) = 7 5 e 
ifei(p) 



po + s a — , p°-4 k+ 

P "'" P 



■A_ 

E+ p 



Po ~ e 



-A+ 

£7 p 



7% 



7% 



f 



-A: 



1 



Pa - ( P 
1 



-A: 



Pa + e? 
1 



Po - ( P 



P0 + (p 



A + A F 



A-A+ 



7 5 £ 



(A-)*A- (A+)*A+ 



£-1 



V A -^A 

with = pg — (e^) 2 — |A ± | 2 , and the three color-flavor matrices are defined as follows: 

= (S ab -S a3 S b3 )5 1 ^ 

(i2)%=s a3 s b3 s^, 

Notice that, when using this quark propagator in loop calculations, one should make the substitutions 



Pa±e p -*p ±e p T «e sign(e p ), 
Po ± £p -> Po ± £p T «e, 



(4) 

(5a) 
(5b) 
(5c) 
(5d) 



(6) 
(7) 
(8) 



(9) 
(10) 

(11) 



in the denominators, and take the limit of vanishing e at the end. This is important for preserving the causality of 
the theory. 
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III. DESCRIPTION OF THE FORMALISM 



Let us start from a simple observation. If the model under consideration had real NG bosons in the spectrum, 
its effective potential as a function of the order parameter AA would have a degenerate manifold of minima. The 
dimension of such a manifold would be equal to the number of the NG bosons. We know, however, that no global 
symmetries are broken in the S2C phase and, therefore, the potential should have a single non-degenerate global 
minimum. The existence of the pseudo-NG bosons means, however, that the potential is nearly degenerate along 
selected directions. The curvature along these directions defines the masses of the pseudo-NG bosons. In the limit of 
the zero curvature, masses go to zero, as it should be for the NG bosons. 

In order to select the directions in the color-flavor space that correspond to the five pseudo-NG bosons introduced 
in Ref . , we should recall their definition. These pseudo-NG bosons correspond to "breaking" of the approximate 
axial color symmetry, given by the following transformations of the quark fields: 

■^d^UV+^d + U^V-^d, (12a) 

4> D 4> D V-U^ + 4> D v+u, (12b) 

U*V-iP% + U T V+^, (12c) 
&+U T + ^V-U*. (12d) 

Of course, this is not an exact symmetry of the model. For example, the kinetic term of the Lagrangian density 
transforms as follows: 

$d (%9 + m° + 4) V>i3 -» i'o (%9 + M7° + P+ujrf + V-rfju) fa, (13) 

and no transformation of the vector field could promote this transformation to a symmetry. 

The axial color transformation, as defined above, allows us to explicitly extract the phase factors of the gap that 
correspond to the nearly degenerate directions of interest, 

A — > V+U^ AU* + V-UAU T , U = exp (iuj A T A ) . (14) 

One could consider ui A as the dynamical fields of the pseudo-NG bosons (which, up to a factor of the decay constant, 
are related to the canonical fields). Such a substitution leads to the following changes of the components of the quark 
propagator: 

iiii -» Rufa) = V+U^RuU + V-URnU\ (15a) 
i? 22 -> R22{u>) = V-U T R 22 U* + V+U*R 22 U T , (15b) 
Ra -> i2i2(w) = V+U^R 12 U* + V-UR 12 U T , (15c) 
R-2i -»i2ai(w) =V-U T R 21 U + V+U*R 21 U\ (15d) 

assuming that the fields lu a are constant in space-time. 

In order to construct the effective potential, we use the C JT formalism jl2j . The corresponding general expression 
reads 

— |jTr \to.G(p)S- x <p) - S- 1 (p)G( P ) + l] + V 2 [G], (16) 

where V2[G] represents the two-particle irreducible (with respect to quark lines) contributions (we will discuss this 
point below). There are, in general, an infinite number of diagrams in V^G]. In our analysis, we leave only a few 
leading order diagrams, graphically shown in Fig. ^. 

Before proceeding to the actual calculation, let us try to understand which type of diagrams could produce a 
non-trivial dependence of the potential on the pseudo-NG fields u> . To this end, we have to recall the origin of the 
pseudo-NG bosons under consideration. In particular, it is crucial that their appearance is related to the breaking 
of the approximate axial color symmetry. It is clear, then, that the dependence of the effective potential on the 
pseudo-NG boson fields results from some mixing between the left-handed and the right-handed sectors of the theory. 
Since there is no any left-right mixing in the diagrams containing a single quark loop [diagrams (a), (6) and (c) in 
Fig. |IJ, the first three diagrams in Fig. |l| turn out to be irrelevant for our calculation. The corresponding contributions 
to the effective potential are free of any dependence on the constant uj a fields. 

Now, we move over to more complicated diagrams (d) and (e) in Fig. |l|. Notice, that these last two diagrams are 
not two-particle irreducible with respect to gluon lines. This is consistent with the fact that we consider the CJT 
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action as a functional of only the quark propagator [for a discussion of this point see Ref. frLql ]. Both diagrams could 
potentially produce non-trivial mass corrections for pseudo-NG bosons. In diagrams (d) and (e), a non-trivial mixing 
of the left- and right-handed quark sectors is possible because there are two separate quark loops. 

A simple calculation shows that the diagram in Fig. [M gives no correction to the effective potential. Thus, the 
leading order corrections come from the diagram in Fig.Tlle. The details of our calculation are presented in the next 
section. 



IV. LEADING ORDER CALCULATIONS 



As we argued in the preceding section, the leading order corrections to the masses of the pseudo-NG bosons come 
from the diagram in Fig. |l|e. In this section, we give the details of the calculation and derive an approximate analytical 
result for the masses. 

The analytical expression for the vacuum diagram in Fig. |l|e reads 

f d^nd^kd^n 

^(Fig. le) = -2 t n 2 a 2 s J P * g Tr [T A »G{p)T B *G{p - q)] Tr [T Al/ G(k - q)T BX G(k)} D fiU (q)D KX (q), (17) 



where the vertex is 



T 



-(T A ) 

The expression in Eq. ( |l7j ) contains two factors of the following type: 

2inas J W? Tl [ rA " G (p) rBKG (p i) 



(18) 



(19) 



In order to extract the dependence of this quantity on the psedo-NG boson fields lu a , we perform the substitutions 
of all component functions, given in Eq. (|l5|). At the end, we expand the result in powers of u) A , keeping the terms 
up to the second order. Thus, we arrive at the following result: 



d 4 p 



tr 



V+rf A Ri2{p)l R (f*) R 2 i( P -q)+ V-rT A Ri2( P h K {T*) T R 2l (p - q) 



V+r{T A ) R2i(ph K T*Ri2(p-q)+P-r[f A ) R2i{p)l K fZR 12 {p-q) 
V + ^f A R 11 {jp) 1 K f B R 11 {p -q)- V-^T^Rniph^Rnip - q) 

V+r (T A ) T R22(ph K (T*) T R22(p - q) - V-r (f A ) T R 2 2( P h K (f B ) T R 2 2( P - q) 



n 



AB,fiK 



X, ,Y fXAC fYBD^jCD.fiK 



(q)+^^r^f 



-U X U Y f XAD fYDC n CB^ (q) + f XBD f YDC u AC^ {q y 



(20) 



Here we used the shorthand notation, 



T A = UT A U f ~T A -lu b f BAC T c + -u b lu c f BAD f CDE T E 

fA = rfj.Au ^ T A +UJ B jBACjC + \^B fBAD jCDEjE 



(21a) 
(21b) 



where f BAD are the structure constants of SU(3). Also, we introduced the one loop polarization tensor in the color 
superconducting phase Q , 



^T A R l2 (ph K {T B ) T R 2 i (p - q) + r {T A ) T R21 ( P h K T B R 12 (p - q) 



7^r A J R 11 (p) 7 K r B i? 11 (p -q)-^ {T A ) T R 2 2( P h K {T B ) T R 22 { P - q) 



(22) 
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By performing the traces over the color and flavor indices, we arrive at the following expression for the polarization 
tensor (see Appendices |X] and for details): 

U AB ^(q)\ AB=123 ^S AB U^( q ), (23a) 
U AB ^(q)\ A B=i 5 6 7 = <5 As nf (q) + i [8 A H m - 6 A5 S Bi + S A6 5 m - 5 A7 6 m ) il% K (q), (23b) 
n 88 >^(g) = ITf (<?). (23c) 

After substituting the expansion ( ppf ) in Eq. ([l7|), we arrive at 

d 4 q 



V(Fig. le) = ij AS-U AB ^iq)Il AB ^(q)D^(q)D KX (q) 
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x | [nf ( g ) - nr (<?)] [nf ( g ) - nf ( q )j + 2nf ( g )nf (<?) 

+ [nf ( g ) - nr (?)] [nf («) - (<?)] | + o [(co A f] . (24) 

It is noticeable that the right hand side of the last expression is independent of the u 8 field, related to the singlet 
pseudo-NG boson. This means that the diagram in Fig. ^e does not give any non-trivial contribution to the value of 
the corresponding mass. Only the ( ant i-) doublet pseudo-NG bosons get a non-zero mass. 

In order to understand this point, it is instructive to consider the "ideal" case, assuming that the axial color 
symmetry generated by 7 5 T 8 is a true (rather than approximate) symmetry of dense QCD. Then, the singlet pseudo- 
NG boson would be related to the breaking of the restricted axial symmetry, defined by the following transformation: 

ip — > exp (iojj 5 ) liip + (25a) 

ijj -> $2i exp (iuj-f 5 ) + $1 2 , (25b) 

^ c -> exp (iu-f 5 ) 1^° + l 2 ip c , (25c) 

jfc ^ ^C Xl exp (iuj-f 5 ) + ^ c l 2 . (25d) 

acting on the first two color quarks. The other symmetry [[^(l)], acting on the third color quarks would remain 
unbroken. This is a simple consequence of the fact that the third color quarks do not participate in the color 
condensation. Now, the axial color transformation generated by 7 5 T 8 could be thought of as the ordinary axial 
transformation f7s(l) accompanied by the unbroken J7s(l). Therefore, it appears to be equivalent to say that cither 
the restricted or the ordinary axial symmetry is spontaneously broken in this ideal limit. In reality, both the restricted 
and the ordinary axial symmetries are explicitly broken. However, while the former is broken perturbatively, the latter 
is broken by much smaller non-perturbative (instanton-like) effects. Because of that, the singlet pseudo-NG boson 
is connected with the ordinary axial transformation U§(1) and its mass is zero in any order of the expansion in the 
coupling constant. 

Non-perturbative analysis reveals that the singlet pseudo-NG boson 77 has a non-zero mass. The value of the mass 
was estimated in Ref. Jf5| . In our notation, it reads 

M^^|A-| 2 expf-^V O^IO 5 . (26) 

The parametric dependence of this mass of the 77 pseudo-NG singlet clearly indicates the non-perturbative nature 
of the underlying dynamics. To derive this estimate, one needs to consider the instanton contribution (notice the 
characteristic exponential factor in the expression above) to the vacuum energy in the color superconducting phase. 

Before calculating the masses of the other pseudo-NG bosons, let us turn to the question of the reliability (and, 
in particular, gauge invariance) of the loop expansion for the CJT effective potential in the present problem. From 
Eq. ( pif ) we see that the vacuum energy is expressed through the one-loop polarization tensor. As is well known, 
the Ward (Slavnov- Taylor) identities (or, equivalently, the BRST transformations) imply that in both Abelian and 
non-Abelian gauge field theories, the polarization tensor is transverse in covariant gauges. We show in Appendix |^, 
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howe ver, that while the one-loop polarization tensor connected with the unbroken SU(2) C is indeed transverse [see 
Eqs. (C3|), ( ]C5| ) and (|Cq)], the one-loop polarization tensors connected with the five broken generators are not [see 
Eqs. (C4), ( Icif ) and ( jp% in Appendix p. 

What is the reason for the violation of the Ward identities in the one-loop approximation? The answer is that it is 
connected with the Meissner-Higgs effect. As was already emphasized in our paper [^) and in the earlier papers Jl6| , p7t , 
this effect implies the presence of unphysical composites having quantum numbers of the (would be) NG bosons in 
any non-unitary gauge, including of course all covariant ones. Although in the unitary gauge these composites are 
"eaten" by the five massive gluons, in other gauges they are crucial for preserving both the unitarity and the Ward 
identities (i.e. gauge invariance) .[] In particular, they lead to an important pole correction in the quark-gluon vertex 
function ||. 

Actually, the traces of this problem are known to appear in other studies in dense QCD. For example, it shows up 
even in the analysis of the Schwinger-Dyson equation in the hard dense loop (HDL) improved ladder approximation 
[^,^8|. In this case, the gauge dependence of the gap appears only through a dependence in the preexponential factor 
of the gap. Formally, it is the next-to-next-to-leading order correction. Usually, therefore, this gauge dependence of 
the superconducting gap is ignored without even studying its real origin. While such an attitude is quite harmless in 
the case of the Schwinger-Dyson equation, it turns out to be crucial in the present problem of calculating the masses 
of the pseudo-NG bosons. Indeed, as is shown in Appendix the gauge dependence of their masses, caused by the 
non-transversality of the polarization tensor, is strong. The reason is that while in the Schwinger-Dyson equation the 
dominant region is that with momenta much larger than the fermion gap |A~|, in this problem the infrared region 
with momenta less than |A _ | dominates. 

This implies that a consistent approximation for calculating the CJT potential has to include those (would-be) NG 
composites. This in turn implies that one should modify the ordinary loop expansion of the CJT potential. Since 
here we are interested in the infrared dynamics, one can neglect the composite structure of these bosons. This leads 
to a non-linear realization of color symmetry breaking. In particular, one has to calculate the polarization tensor in 
this framework. 

In Appendix |^, we consider this problem. As is shown there, including the (would be) NG bosons indeed restores 
the transversality of the polarization tensor. The remarkable thing is that the dominant contribution to the masses of 
the pseudo-NG bosons c ome s only from the part connected with the magnetic gluons [singled out by the projection 
operator see Eq. (|Cl|)], and this part, unlike the other contributions, remains unchanged by including the 

contribution of the would be NG composites. This implies that for the calculation of the masses of pseudo-NG bosons 
in the Landau gauge one can use the initial, unmodified framework for the CJT effective potential (at least in this 
approximation). Indeed, the dangerous longitudinal terms in the polarization operator does not contribute in this 
gauge and only the contribution of the magnetic gluons matters. Because of this observation, henceforth we will use 
this gauge. 

To get a rough estimate of the m ass o f d ouble ts, on e cou ld use the following order of m agnitude asymptotes for the 
polarization tensor [see Eqs. (p8|), ( |BI4D , ( |B20[ ) and ( |B24| ) in Appendix ||] and Eq. ( |Cf| ) in Appendix |cj: 



ni 1 ^) 



n?>( 9 ) 



a at i, 2 \A-\ 2 
q s ^ 2 |A-| 2 

q 



In 
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for |g 4 U«2|A-|; 

for g,|A-| < M; 

for | A- 1 < |g 4 | < q- 

for |g 4 | < | A" | < q. 



(27) 



By making use of these asymptotes along with the HDL expression for the gluon propagator [see, for example, Ref. 0], 
we check that the dominant contribution to the quadratic term in Eq. ( p4[ ) in the Landau gauge comes from the region 
of momenta where |A| < q < [i and < | £74 1 < |A|. Therefore, by taking the explicit expressions for the polarization 
tensor given in Appendix |c] into account, we derive 



M 



1 d 2 V I70tt 



F 2 du 2 
340tt 



-ot s 



3fi 2 



q 2 dq 



|A| 



|A| 



ln IAI ~ Cm ^* s 



dq 4 



{q 3 + 
Cm 



5™d 



9*1) 

857T 2 



a s ^ 2 \A\ 



V2^ « 7 x 10 2 



(28) 



1 Note that, because of the composite (diquark) nature of the order parameter in color superconductivity, it does not seem to 
be straightforward to define and to use the unitary gauge in this case. 
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where the definition of the Debye mass, m 2 D = (2/7r)a s /U 2 , the decay constant F 2 = fi 2 /8ir 2 , as well as the relation 
between the gap and the chemical potential, ln(ju/| A|) ~ 3(/7r/2) 3 / 2 aJ , were used. 

V. PSEUDO-NG (ANTI-) DOUBLETS VERSUS CONFINEMENT 

In the color superconducting phase in the model with two quark flavors, the SU(2) C subgroup remains unbroken. 
Since the only massless quarks of the third color do not interact with gluons of SU (2) c , the corresponding dynamics 
of gluons decouples in the far infrared region, p -C | A~ |. The low energy effective action of the SU (2) c gluondynamics 



was derived in Ref. 19 1. Of special interest for us here is the observation of Rcf. [|19| that the unbroken SU(2) C is 
confined at sufficiently low energies. The corresponding scale was also estimated. It is given by the following simple 
expression: 

A^ CD ^|A~|exp[-Coa 2 exp(C/^7)], C = 3(|) 3/2 , (29) 

(there is some uncertainty in determining Co — 10~ 3 ). 

Now, let us discuss how such confinement could affect the properties of the pseudo-NG bosons. In particular, this 
concerns the doublet and antidoublet states which are colored under the leftover SU (2) c subgroup. Because of the 
confinement, these pseudo-NG states cannot exist as free particles. Instead, they should form some colorless bound 
states. 

The colorless bound states should be somewhat similar to charmonium in ordinary QCD. In fact, the situation with 
the doublct-antidoublet states in the color superconducting phase is much simpler. This is because the new confinement 
scale Aq CD is extremely small as compared to the masses of the (anti-) doublets themselves. This already suggests 
that the binding dynamics is essentially perturbative [this could also be checked posterior, see Eq. (|32]) below]. The 
bound colorless states would be similar to the positronium in QED pp|-p2[|. The only difference is due to the color 
superconducting medium effects. The Coulomb potential between two static charges at a distance r from each other 
is given by |L9[ 

47TQ; I I 

W\ <<r<< ^ m 



where 



(31) 



9tt|A-| 2 

is the dielectric constant [ fl9| [see also our derivation of Eq. (C3) in Appendix 0. By analogy with the positronium 
(notice that the constituent doublets are spinless in the problem at hand, but such a difference affects only the fine 
structure of the spectrum), we get the following estimate for the binding energies of the colorless states built of the 
doublet- antidoublet pairs: 

Ma 2 „ a s - 39/4 |A"l ( Z^'A 
E n ^—rr%- c B— J [ ex P ~ j— )' n=l,2,..., (32) 

with C E ~ 1.8 x 10 11 . Notice that this binding energy is parametrically much larger than the confinement scale 
A.'q CD in Eq. ( ^9| ) when a s — > 0. This means that it is the (perturbative) Coulomb contribution to the potential 
that is mostly responsible for the paring dynamics of the (anti-) doublets into colorless hadrons. The linear confining 
(non-perturbative) part of the potential is of minor importance at least for the low lying states. 

We recall that in the S2C phase, there exists a conserved (modified) baryon charge. While the value of this charge 
equals zero for the quarks of the first two colors, it is +1 for the third color (massless) quarks. Since the pseudo-NG 
(anti-) doublets are composed of one massive (anti-) quark and one massless (anti-) quark, they are scalar baryons, 
carrying the baryon charge ±1. The colorless hadrons composed of them (discussed above) can be called mesons. 

VI. CONCLUSION 

In this paper, we derived analytical estimates for the masses of the pseudo-Nambu-Goldstone bosons in the color 
superconducting phase of dense QCD with two light flavors. In agreement with our previous hypothesis, the mass 
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values are small compared to the value of the superconducting gap at sufficiently large quark densities. Analytically, 
the expression for the masses of doublet and antidoublet pseudo-NG bosons are suppressed by a power of the coupling 
constant. The mass of the singlet is exponentially suppressed. The mechanism for mass generation of the singlet is 
purely non-perturbative m| . 

As was shown in fll9| |, the unbroken SU(2) C subgroup is subject to confinement. This fact has an important effect 
on the properties of the doublet and antidoublet pseudo-NG bosons which are colored with respect to SU(2) C . In 
particular, we argue that the doublets and antidoublets should form charmonium-likc colorless bound states. Moreover, 
since the confinement scale of SU(2) C subgroup is extremely small compared to the value of the masses of doublet 
and antidoublet, the pairing dynamics of the colorless mesons is dominated by perturbative Coulomb-like forces. This 
allows us to estimate the corresponding binding energy. 
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APPENDIX A: GENERAL FORM OF THE POLARIZATION TENSOR 



The general form of the polarization tensor in the color superconducting phase was given in |Hj . For completeness 
of the presentation, we also rederive it here 



n^(s)L, B=1 , 2 , 3 - ^ITO = 2i5 AB ,a s J 



x tr D 



\E+(p) E-^p)) 1 \El{p-q) E+(p-q) 

+ y / (A~)*A. p (A + )*A+ \ / A + A~_ q A-A+ \ 
V E A (p) E+(p) J 1 \E+(p-q) E A (p-q) 



+ 7^7° 



A ' 

Po +e: 



-A: 



Po - e 



E A (p) p E+(p) 
-A 



^A+ 



( Yl o ( PQ-1Q + * P - q 



E A (p-q) 



A 



p-q 



E+{p)^ ' E A (p) 



Po- qo - £p- q A+ 



Et(p-q) 



p-q 



p -q -e p _ 
E A (p-q) P 



(Al) 



U AB ^{q)\ A B=4 5fi 7 = S AB U^(q) + i (S A4 S B5 - S A5 S B4 + 6 A6 5 BT - S A7 S m ) fif (q) = nra s 6 



. IB 



d 4 p 



(2nY 



+ 7^7° 



tr_o 



1 



7 M 7° 



-A: 



Po - e P Po + e-p 



P° + e P a- , PQ-S A + 
E A (p) ? + E+(p) p 



7^7° 



-A: 



-A+ 



P-Q i + P—Q 

po- qo- e p - q Po - qo + £p- q 



+ 7 M 7 U 



( P0 + 4 A - j_ P°- £ P A +^ ^„,0 



E+(p) p E A (p) 



1 



-A; 



1 



-A+ 



Po - £p Po + e P 



p 7 7 



Po - qo + 


e p-q 


E A (P- 


q) 






Pa - qo - 


e+ 

P-Q 


Po~qo + 


f+ 
P-Q 


Ei(p- 


q) 



A; 



A: 



El(p 


-q) 


1 




po - qo 




Po - qo 


_ e p-q 


E A (P 


-q) 



^Aj 



a; 



A _ ^ — ^ PZ1 A + 

p-q ' jp—f m „\ p-q 



(S AA 5 B5 _ S A5 S B4 + S A6 S B7 _ ^7^6) 
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d 4 p 



„„„ ( PO + £p A - | PO ~ 4 



Tl 



1 



i 



A+ ) 7 " 7 ° 



1 



-A„_„ + 



-A+ 



+ —^ A p ^7 



Po - e P Po + e P 



7 M 7° 



-A: 



-A+ 



+ p i — p 
PO - £p Po + e P 



7 7" 



PQ — (JQ + 


e 

p-q 


E A (P- 


q) 






Pq-Qo- 


e+ 
p-q 


po-qo + 


e+ 
p-q 


eKp- 


q) 



Po-qo- £ p -q " - Po - qo + e 
Po-go-e+_ g + 

P-9 " t " 77+^ „\ P-9 



A P + - 9 



P-9 



eKp- 


q) 


i 




Po - qo + 


e p-q 


Po - qo - 


e p-q 


Ea(p~ 


q) 



K- q 



(A2) 



x tv D 



Y f A+A p + r A a ~)* a p-. x ( a+ )* a pV 



(A-)*A- (A+TA+ ^ ^ f A+A p _^ 
££(p) 25+ (p) 



A- A 



g-g 



7"7 



E A ( P ) p+ E+ip) v) 1 ' \ E A (p-q) " ' 



££(p-g) E A {p-q) i 
Po - qo + e~ 



<-A+ 

Etip-q) 



u o ( p° + 4 a - , pq- £ p A +A „ o ( Po — qo + e p-q A - - '/n 



A" + — - i— 1A 



- 27^7° 



1 



1 



Po - e P Po + £p 



_A " + C + A p hi' L „ ^- •■/• / ,. 



-A„_„ + ■ 



-A+ 



27^7° f ^-^A p " + — ^-A+) 7 " 7 ° 
VPO - £p Po + £p / 



v Po - 9o - 
1 



cp Tt p / \P0 — qo — e p-q 

Here we made use of the the following results for color-flavor traces: 



— K- q + 



Po - q + e p _ q 
1 

po - qo + £p- q 



A P% 



tr c/ (l.T^T 3 ) 
tr c/ (I 2 T A I 2 T B ) 
tr c/ (liT A I 2 T B ) 



S AB A,B = 1,2,3; 
A, B = 4,5,6,7; 
| A,B = 8: 

i4,B = 1,...,7; 
| .4,5 = 8; 



1 xAB 



q -y 2 3' 

5 AB + i (^4^5 _ 5 A5 5 B4 + _ ^7^56) ^^5 = 4^ 5^ 6^ 7; 



tr, 



c/ 



iT A i (T B ) T 





-<5 AB ^,5 = 1,2,3; 

A, B = 4,5,6,7; 

1 A,B = 8. 



A,B = 8: 



(A3) 



(A4) 
(A5) 
(A6) 

(A7) 



APPENDIX B: CALCULATION OF INTEGRALS 



In this Appendix, we calculate different type of integrals that appear in the expression for the polarization tensor 
(see the previous Appendix). For our purposes it is sufficient to consider the gluon momenta much less than /x, 
neglecting all corrections of order g 2 /p 2 . In this approximation, the Dirac traces, containing the on-shell projection 
operators, read 



9 



tr 



tr 



7 , A ( ± y A ( ± _) 



tr 
tr 



IpT J ' 

IpI 



2[g ,o g ,o_lf_ 



IpI 

IIP 



(Bl) 
(B2) 
(B3) 
(B4) 



Also notice that in order to perform the integrations properly, one should make the following replacements in the 
denominators of propagators: 



E% - E% + ie, 



Po±e p -> po ± e p T «e sign(e p ), 
Po ± £p -> Po ± £p T ie. 



(B5) 
(B6) 
(B7) 



Form the definition in the previous Appendix, we see that the general structure of the polarization tensor reads 

nn?) = 4^ s (J^(q) + I^{q)) , (B8a) 



nf (g)=47ra 5 /n<7), 

nrc?) = (^(?) - (9)) + ^^sw?)- 



[Here we took into account the fact that the off-diagonal term in II j4B ' AU/ (g) 
to calculate the corresponding integral in Eq. ( A2).] Let us start with the calculation of the first type of the integrals, 



A, 3=4,5,6,7 



(B8b) 
(B8c) 

is zero. To see this, one has 



HDL 



(<?)=i 



tr 



= -2 



d 4 p 
(2rr) 4 

d A p 
(2^ 

d 3 p 



D 



tr_o 



7 M 7° 



7^7° 



-A; 



Po - e P Po + £p 



-A; 



-a; 



p— i i + p—q 

Po-Qo- e D - Po - 9o + C-o , 



~^ A p + ^- A p ) 7 l/ 7° 
Po - £p Po + e P 



-A; 



-a: 



+ p—q i — p—q 

Pa- qa- £p- g Po - ?o + e p _ 9 



(2rr)3 
d 3 p 



0(-Cp)fl(^L 9 ) 0(e;)0(-e;_ g ) 



P-9 



g + fip-q - £p + go + ie 



(2tt)3 



IpT 



p p 



?0g f </u 90 / ((/ „,^ 



92 g 2 



,i/0 9 



5" 1-3 

9 g 







(B9) 



where 



'(*) 





1 + » 






- 1 




1 -X 





• 7r l 

'2 1 



(BIO) 



It is easy to check that q^I l ^ l/ DL (q) = 0, as it should be. Also notice that in Euclcdian space, x — iq^/q = iy, and 



Q{x) Q (y) = — 1 + y arctan - 

y 



■l + fy, for y<l; 
-3^, fory»l. 



(Bll) 



Now, let us consider 
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2i 



d 4 P 

d A p 
d 4 p 

2 r 1 



tr D 



tr D 



/A+A- A-A+A f (A-yA-_ q < (A+)*A+_; 



fi^N |A~| 2 



p 2 y E A (p)E A (p-q) 
i I A- 1 2 [(1 + C 2 ).g^V° + (1 - C 2 )3^ + (1 - 3^ 2 )^ 



\A-\ 2 + x{i-x) {eq 2 -ql)-ie 

After switching to the Eucledian momenta (qo = iqi), we get 



(B12) 



i a- i 2 (i + a^v + (i - + (i - 



\A-\i+x(l-x) (qi + eq 2 ) 



i ia-| 2 [(i+^w^a-aff^ + a-sa^ 



V^+e 2 gV94+C 2 9 2 + 4|A-| 2 



x In 



v /g 2 + ^V + 4|A-| 2 + y g4 2 + e 2 g 2 



Vii + eq 2 + 4\A-\ 2 -Vii + eq 2 

In different limits, we obtain the following behavior: 



(B13) 



3tP 



■TO = 



2M 2 |A | 2 , g| /„m0„^0 , 1 „/«A 



for |g 4 | l? «2|A- 
for g,|A-| < M; 



3tt5^ ^ |A- 

-S'^Afj""^^/^^, for \A~\«\ qi \«q; 



iK\qi\q iix JA 1 
2 |A- 
8<? 



for |<j 4 | < | A" | < g, 



where (in Minkowski space) 



(B14) 



W (2) (g) = 



90tt 2 |A- 



(B15) 



As one could see, the imaginary part of J A '(q) in Eq. (B12) appears only for qo > 2|A |. Its explicit form is given 
in terms of the elliptic integral of the first, F(ip, z), and second kind, E(tp, z), and the hypergeometric function, 



lm[J^(q)} 



M 2 |A" 



-g 2 -4|A- 



27T ^q2 ^ q 2 _ 4 | A -|2 

q z 



F arcsin < 



g 2 



M 2 |A-| 2 0(4|A-| 2 + g 2 - g 2 )0( g 2 -4|A 



F ^arcsin 

- 12\ 



g 2 



q 2 o q 2 o~MA- 



2tt 

Trfag -4|A-| 2 ) 
4<? 2 



g 2 Vgf 



q 2 

1 3 g 2 -4|A-| 2 
2' 2'"' g 2 



A' 



E ^arcsin 
q$ - 4|A" 



9o 2 '<7o 2 -4|A-| 2 

2 



g 2 



% 2 9o 2 -4|A- 



g 2 



(B16) 



Similarly, let us calculate the following quantity: 



d 4 p 

d 4 p 

(27T) 4 



tl'D 



E A (p)' Ap ' 



'-A 



Po-go-£p- g . + 



7^7 



o ^ PQ + 6+ A _ , po - e p ^ ( Po-qo + e p _ q A _ , Po-go 



■A" + — -£-A + 



^A(P-ff) 



-A+ 
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3tt 2 



d 4 P 

(2tt)4 



IpI 



P^A (Po + £p)(Po~gQ + £ P-9) 



IpI 



E A (p)E A (p-q) 



- 3?r2 (5 9 9 )+2*J o dxj (27r)4 + \g + ^ 

[po + £p + x(q + £q)] [po + e~-(l- x)(g + 



- (ep ) 2 - x(l - x)(eq 2 - ql) - |A" | 2 + ie 



2 ' 



(B17) 



So, we derive 

TO 



\A-\ 2 + 2x(l-x)(q qt + q*e) 



+ x(l-x)(<r?-<$)-ie 



1(3 - e W° - l(i - £ 2 )<r - l(i - 3e 2 )^f + e f + 



3tt 2 



(<T - ^V°) + ^ J dxj dC (3 - e 2 )^°^° - (1 - (1 - 3C 2 )^| 

|A-| 2 + 2x(l - x)g 2 ^ 2 



lA-p+aKl-aOGr^-ggJ-fe 



(i 



2 /•! 



Jo 



+ ^ rfx d£ ^ u - + <r- 



9 



a;(l - a;)^^ 2 



g/ |A-| 2 + .T(l-.T)(g 2 e 2 -^)-ze' 



(B18) 



After switching to Eucledian momenta (g = iqA), we could also perform the integration over x, and get the following 
result: 



TO =I%DL{<l) + 7^j\ 

\A-\Hqj-eq 2 ) 



(3 - e w° - (i - £ v a - sa 



Vg + eg 2 + 4|A~| 2 + y/g 2 + gy 



■In 



2^ (^°? + ^°f)|A-|W 

—5 / 



■In 



y/gl+eV + 4|A-| 2 + y/gl+eV 
^+CV+4|A-|2(g2 + e2(? 2 )3 /2 " v /g 2 +e 2 g 2 + 4|A-| 2 - + £V ■ 

The following asymptotes are valid 



(B19) 



l£ v (q)-lHDL(q)=< 



^ In ^ (y*** - 5 - + 2^f^) 



for |g 4 |,g«2|A-|: 
for g,|A-| < |g 4 |; 



" 47r| g4 |g m |A-| = 



where (in Minkowski space) 51^ (q) reads 



i/0 _ a ixv _ <r q 



(B20) 



« (2) (g) = 



90tt 2 |A-| 2 

Finally, we calculate the last type of integrals 
d 4 p 



^V°(3g 2 + Wq 2 ) - -g^(q 2 + hq 2 ) + <t<? 



+ 



for \q 4 \ < |A"| < q. 



fi 2 q {q^u v + u^q v ) 
18tt 2 |A-| 2 ' 



(B21) 



TO = 



(2tt)4 

d A p 
(2tt)4 



trzj 



tr_c 



7^7° 



7^7 



A" + 



Pa - £p 




7^7° 




1 






^Po 


-Qo- 


e p-q 


Po - £p 


A+" 


7^7° 




1 








- 9o - 


e+ 
p-q 



a: 



-a: 



Po - <7o + e p _ q 



3tt 2 



(5"" - 9^9 V ° 



if 



PO - 50 + tp-q 

(PO +£p )(P0 - 90 + Ep-g) 



A+ 

p-q 



\P\J [p 2 -(ep) 2 -\A-\^}E A (p- q ) 
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[Po + e~ + x{q + £g)] [p + e~ - (1 - x)(q + £g)] 
[pi - (e p -)2 - x (l - z)(£V - g 2 ) - (1 - x)|A-|2 + fe] 



g*° + ^- 



3tt 2 
, A* 1 



4tt 2 



(3 - a^V - (1 - e)9^ - (1 - 3£ 2 )^£ (I A- 1 2 + 2xg 2 ^ 2 ) 



|A-| 2 +:r(g 2 £ 2 -g 2 )-z£ 







9 



xqqo£ 2 



qj \^+x(q^e-ql)-^e- 



(B22) 



After switching to Eucledian momenta (go — *94)i w ^ could also perform the integration over x, and get the following 
result: 



(3 - £VV° - (i - £ V - (i - 3£ 2 )^T 



M 2 «<?4 



71" 9 JO 



/ d ag> iv — + g' j0 — 



2„2 



(q 2 + eV) 2 



94 + r g 



|A- 



|A-| 2 (gi-eV) ln LA-| 2 + (? 2 + ev 



{ql+iW 



|A- 



(B23) 



The following asymptotes are valid 



-iHDL^ + ^ig^g^ + h* 1 ") 



I^{q)-I HDL {q) = { 



3tt 2 



w (2) (g), 



for |g 4 |,«<2|A-|; 



(^"V " <T + 2 ^V+fy° ) , for g, |A~| « |g 4 |; 



87r|g 4 



47rg 



1 (w°g 



f?4 



I/O _ M" _ 



+ % 



for I A- 1 < |q 4 | < g; 

for |g 4 | < | A" | < g. 



(B24) 



where (in Minkowski space) 



SI (2 \q) 



30tt 2 |A- 



p 2 q (g^°+g^V) 
5tt 2 |A-| 2 



(B25) 



APPENDIX C: GAUGE INVARIANCE AND POLARIZATION TENSOR 



By making use of the definitions in Eq. (|BS| ) as well as the general structure of the tensors J^(q), 1^ '(g) and 
l£ (g), given in the preceding Appendix, it is easy to show that all three polarization tensors Tl^ 1 " '(g) (« = 1, 4, 8) take 
the following general form: 

nn?) = nf ) (g)o( J )^(g) + n 4 (2) (g)o( 2 )^(g) + nf \q)o^ v (q) + nf > (g)OW^(g). (CI) 

In this representation, we use the following set of four tensors p3] : 



C$(g) = g Mi / + *L , 


(C2a) 


(2) (o) -u •• ^ ^ 


(C2b) 




(C2c) 


^(9) = oiV| + |^- 


(C2d) 



The first three of them are the same projectors of the magnetic, electric and unphysical (longitudinal in a 3+1 
dimensional sense) modes of gluons which were used in Ref . |7(] . In addition, here we also introduced the intervening 
operator 0${q), mixing the electric and unphysical modes p3| . Note that u M = (1,0,0,0) and q^ = g M — (u ■ q)u^. 
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By making use of the representations for different types of the polarization tensors in Appendix [b|, we derive the 
following explicit infrared (|go|,g |A _ |) asymptotes: 



nrc?) - 



9tt|A-| 2 
9tt|A-| 



[rfVV + g 2 ) q 2 ^ + q a (g*V° + S^V)] 



q 2 O^(q) + (q*-q*)O^(q) 



3 „it.u / \ 2,cy s fj, 



nf (g) = ^nf (g) = 



3tt 



(1) ""(g) + 



<Zo-3g 2 ^ f2W _, , 3g, 2 - 



<?o ~ I 2 



, (2)^( g ) + 



2 2 

95 - <r 



0(3)/"' (g) 



(C3) 
(C4) 



Similarly, we derive the asymptotes in other regions of interest: 



■ 8a t4~' ^]Z^[0 (1) (<l) + OW( q )], for g ,|A-|«|g 4 |; 



V.>?{q)-^a s I% DL {q) = { -B^pi^^oWiq) 



|A-| 2 



2wa ^ |A '' 0( 2 )(g), 



for | A" | < |g 4 | < g; 
for |g 4 | < | A" | < g, 



(C5) 



nr(9)-47ra a 7^ M («) 



2a sM 2 |A-| 2 , 

i — 1 2 



2|g 4 |g 

7TOt s fl 2 | A~ 



, L o«(g) + o( 2 >(g)- 

(1) (?)-0 (3) (?) + t0W(g)' 

0(l)(g)_0(2)( g ) + iS4 (4)( ? ) 



30<%) + £oW(g) 



for g,|A"| < \ q4 \; 
for | A" | < |g 4 | < g; 
for |g 4 | < | A- 1 < g, 

(C6) 



and 



nf(«)- 47TQ s 7^ L (g) 



3^ 
2a a M 2 |A~ 

3|?4|g 

7ra s /j 2 | A - 

3g 



^ln 



it 



-ln^L 

m |A-| 2 



0( 3 ) (g) 
OW(g). 



i|iQ( 4 )(g) 



• *2_0W(<i)' 



for g,|A-| < |g 4 |; 
for | A" | < |g 4 | < g; 
for |g 4 | < | A - 1 < q. 



(C7) 



As is easy to check, II^ 1/ (g) is transverse everywhere, while Hj L l '( q ) and IIg 1/ (g) are not. From the structure of the 
polarization tensor H^ u ( q ) in Eq. (C3), we derive the explicit expression for the dielectric constant in the far infrared 
region where the three gluons of the unbroken SU(2) C decouple from the other degrees of freedom. It coincides with 
the result in Ref. which is quoted in Eq. (pi). We also see that the mag netic permeability is equal 1, because 
there is no magnetic contribution g 2 ( - 1 - )A "'(g) on the right hand side of Eq. (C3) [notice that q^O^ 1 ^ 11 (g) is the electric 
type contribution]. Let us also notice that the general expression for II]" y (g) is transverse for all momenta. To see 
this, we derive the following representation: 



U^(q)=^aJ^ DL (q) 



4a s /x 2 |A" 



1 [g?(l ~ e)Q {1) ^{q) + 2£ 2 (gf + g 2 )Q( 2 >Hg)] 
V<74 2 +£V+4|A-|2(g| + £V) 3/2 



x In 



y/q\ + £ 2 g 2 + 4| A^F + y/qj + £V 
v/g 2 + C 2 g 2 +4|A-|2- ^ql + eq 2 ' 



(C8) 



From the asymptotes in Eqs. (B14) and (B2C), one might get the impression that the transversality of n^^g) = 
4ira s [ J^(q) + / MI/ (g)] is not exact. This is just an artifact of using the expansion for either the limit of |g 4 | <C g or 
g <C |g 4 |. To test the condition g Ai n^ 1/ (g) = 0, one has to keep the subleading corrections to such expansions. 

It is straightforward to check that the general expression in Eq. ( ^4|) becomes strongly gauge dependent if the 
one-loop polarization tensors n 4 Ily (g) and H^ u (q) with nonzero longitudinal components are used in the calculation. 
Moreover, the corresponding gauge dependent contribution to the value of the mass comes from the region of small 
momenta, |g 4 |,g <C |A _ |, and it is logarithmically divergent, 



\ 2 a 2 a ii 2 In ■ 



|A-| 



(C9) 



where A is the gauge fixing parameter and em is the infrared cutoff. In deriving this result, we used the asymptotes 
presented in Eq. (04). This result clearly demonstrates that the presence of the longitudinal components in the 
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polarization tensors U^ v (q) and IIg 1/ (g) is not acceptable. In the complete theory these components should be zero, 
because the Ward (Slavnov- Taylor) identities imply that the polarization tensor should be transverse. As we mentioned 
in the main text, it is the contributions of the (would be) NG bosons which should be added to the polarization tensors 
in order to restore the transversality. 

As it was already pointed out in Sec. one has to use the non-linear realization of color symmetry breaking for 
that purpose. The polarization tensor can be found from the corresponding low energy effective action (a similar 
approach for improving the polarization tensor has been recently considered in Ref. 23 1). The effective action for the 
two flavor dense QCD was described in Ref. It is given in terms of the Maurer-Cartan one-form, 



w M = VH(dp - igA^)V, 

where the field V parametrizes the coset space SU(3) C x U(1)b/SU(2) x U(1)b- Its explicit form reads 



V = exp 



A=4 



(CIO) 



(Cll) 



In the last expression, we took into account that one of the broken generators is (/ + /3 rather than T 8 . The 

fields <j) A (A = 4, . . . , 8) describe dynamical (would be) Nambu-Goldstone degrees of freedom which should necessarily 
appear in the low-energy description of the color superconducting phase of dense QCD (unless a unitary gauge is 
used) . The interaction of gluons and would be NG bosons with fermions is described by the following term |24|] : 

L / =^ + |i7°+y i w |t )^. (C12) 

Therefore, after integrating out fermions and high momenta gluons (with q 3> \q±\ |A _ |) which are mostly re- 
sponsible for generating the gap , the one-loop approximation of the effective action is mimicked by the following 
expression: 



S eff = -o Trln 



K.+ A 
A £_ 



--Tr In [/C_/C+] — ^Trln 



l-JCT^A/C^A 



where 



lC+ = + /i7° + VH{$ - ig#)V, 
= ^7° + V T i{$ + ig4 T )V*. 
The quadratic part of the induced effective action, as follows from Eq. ( |C13| ) , reads 

d 4 qd i p 



(C13) 

(C14) 
(C15) 



S 



(2) 



-tr 



Ru (p) f+(q)Rn (p - q) /+(-?) + Rl2 (p) f-(q)R21 (p - q) /+(-<?) 



+ R 21 (p) f+{q)Rx2 (P - q) f-(-q) + R22 (p) f-(q)R22 (p - q) /-(-<?) 
where, by derivation, 



/+(<?)= 7 A 



/-(<?)= 7 A 



gA A (q)T A + l q^ ( t ,A (l) TA +^'- 



'(?) 



V3T 11 



-gA A (q)(T 



A\T 



1% 



1% 



A=4 



(C16) 



(C17) 



(C18) 



(here we expanded V in powers of <p A ). By substituting the last expressions into Eq. (C16), we arrive at the result in 
the following form: 



c(2) 
a eff 



d 4 q 
WT 4 



J2A A (-q)U^(q)A A (q) 



A=l 



A=4 



A A {-q)-- g q^ A {-q) 



AU-q)-^-q^(-q) 



nf(<?) 



A A {q) + -q^ A {q) 
9 

AUq) + —q^(q) 
9 



(C19) 
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where, by derivation, 



U^(q) = - [lL^(q)+8na s I^ DL (q) 



(C20) 



By integrating out the cf) A fields, we arrive at the effective action of the gluon field [compare with the effective action 
in Ref. ||]: 



•J, 



(2) 



d 4 q 



(2tt)4 



A=l 



nf (g) - Uf(q) qit , (qxli^q.y 1 q v ^t (?) 



A=4 



fir (g) - nr'(9)v (^n^) 1 ^nr (?) 



(C21) 



By making use of the properties of the polarization tensor, we check that this quadratic part of the action splits into 
three decoupled pieces with different types of gluons. One of them contains the first three gluon that do not feel the 
Meissner effect. The improved polarization tensors for the five gluons corresponding to broken generators, explicitly 
transverse, and thus, they can be written in the standard form: 



? 9 

<7o -T 



n M (<?)0 (2) ^(?), 



(C22) 



where 



U itt (q) = -UP(q) 



Ri,i(q) = — 2 2 

9o -T 



nf( 9 ) 



[n< 4) (g)] 2 
n| 3) ( 9 ) 



(C23) 
(C24) 



with ll[ x) (q) (i = 4, 8 and x = 1, 2, 3, 4) defined in Eq. (JcTJ) . 

In the infrared region |g 4 |, q <C | A - 1, the improved gluon tensors that correspond to the broken generators read 



3tt 



0^(9)+ t f 9 0^(q) 



(C25) 



Similarly, we get the expressions in other limits: 

2a s ^ 2 |A~ 



3tt(J_ : 



hi 



?4 



hdl(i) 



|A-| 2 



[Oa)(g)+0(2)( 9 )], for g ,|A-|«M; 



^^[0^(q)-0^(q)] 

TO .^|A-| [ (l)(g)-Q(2)(g)] 



for | A" | < |g 4 | « g; 
for |g 4 | < I A- 1 < q, 



(C26) 



and 



n 8,ne«,(9) " 4Tra s l£ v DL (q) 



0, for g,| A- 1 « |g 4 |; 

^!iA_Pln^OW(g), for |A-|«|g 4 |«g; 

m, f ' ( "(g), for | 94 | « I A- 1 « q. 



(C27) 



One should notice that the magnetic components [determined by the terms with the projection operator (q)} of 
the new polarization tensors, are exactly the same as in the one-loop approximation given at the beginning of this 
Appendix. The electric components are different in general, but they appear to be also essentially the same to the 
leading order approximation. Most importantly, however, the longitudinal contributions in II 4 (g) and H^ v (q) are 
gone now. The immediate consequence of this is that the vacuum energy given by Eq. (£4j), as well as the masses of 
the pseudo-NG bosons, become explicitly gauge invariant. 

For completeness of presentation, we also note that the contribution of the modified electric components of the 
polarization tensors to the value of the pseudo-NG boson mass is of order 



1G 



5 el M 2 ~ l - L (C28) 

A 4 

Because of the chemical potential in the denominator, this is clearly a subleading correction as compared to the result 
in Eq. (ph. 
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